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Inertial Range and the Kolmogorov Spectrum of Quantum Turbulence 
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Quantum turbulence is numerically studied by solving the Gross-Pitaevskii equation. Introducing 
both the energy dissipation at small scales and the energy injection at large scales, we succeed in 
obtaining the steady turbulence made by the balance of the injection and the dissipation of the 
energy. The inertial range takes the Kolmogorov spectrum and the energy flux is consistent with 
the energy dissipation rate at small scales. These results confirm the properties of the inertial 
range of quantum turbulence, proposing a prototype of turbulence much simpler than conventional 
classical turbulence. 

PACS numbers: 67.40.Vs, 47.37.+q, 67.40.Hf 



43 

o 

0> 



-a 
c 

o 
o 



> 

On 
O 

o 
in 
o 

-I— > 
c3 



i 

C 

O 

o 



Quantum turbulence (QT) consisting of quantized vor- 
tices has become a great subject not only for investigating 
the issues of superfluidity in the proper field of low tem- 
perature physics but also for studying a prototype much 
simpler than conventional classical turbulence (CT). The 
steady state of fully developed CT of an incompressible 
fluid takes the characteristic statistical law for the energy 
spectrum 0]. The energy, injected into the fluid at some 
large scales in the energy-containing range, is transferred 
in the inertial range to smaller scales without been dis- 
sipated. In the inertial range the energy spectrum takes 
the Kolmogorov law 



E{k) = Ce 2 / 3 £T 5 / 3 . 



(1) 



Here the energy spectrum E(k) is defined as E — 
J dk E(k), where E is the kinetic energy per unit mass 
and k is the wave number from the Fourier transfor- 
mation of the velocity field. The energy transferred to 
the energy-dissipative range is dissipated by the viscos- 
ity with the dissipation rate e of Eq. Q , which is equal to 
the energy flux II in the inertial range. The Kolmogorov 
constant C is a dimensionless parameter of order unity. 

The inertial range is believed to be sustained by the 
self-similar Richardson cascade process, in which large 
eddies are broken up to smaller eddies through recon- 
nections. However, the Richardson cascade has never 
been confirmed clearly in CT, because it is impossible 
to definitely identify each eddy. In contrast, quantized 
vortices are definite and stable topological defects with 
the quantized circulation, consequently QT may give us a 



prototype of turbulence much simpler than CT |2( ■ Some 
experimental 0, 3 or theoretical works 0, @ studied the 
Kolmogorov law for QT, but it was very difficult to con- 
trol the energy dissipation rate e. 

Recently we made the numerical simulation of the 
Gross-Pitaevskii (GP) equation which describes a quan- 
tum fluid at zero temperature |2|. Here we introduced 
a dissipative term that works only in the smaller scales 
than the healing length, in order to remove compressible 
short-wavelength excitations created through vortex re- 
connections which may hinder the cascade process. We 
thus obtained the Kolmogorov spectrum for decaying tur- 
bulence starting from the random configuration of the 
velocity field. In this work, we also inject energy into 
the fluid at large scales by moving the random potential, 
thus obtaining steady turbulence sustained by the bal- 
ance of the injection and dissipation of the energy. The 
energy spectrum of the incompressible kinetic energy still 
takes the Kolmogorov law, and the energy flux n is al- 
most independent of the scales in the inertial range and 
consistent with the energy dissipation rate in the energy- 
dissipative range. These results make clear the properties 
of the inertial range of QT. 



THE MODEL OF THE GP EQUATION FOR 
STEADY TURBULENCE. 



We study numerically the Fourier transformed GP 
equation 
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-^£v(fci,t)*(fc-fci,i). (2) 



Here $(fc,t) is the Fourier transformation of the macro- is given by 
scopic wave function <!>(x,t) = f(x,i) exp[i(f>(x,i)}, g is 

the coupling constant, L the system size. The dissipation j(k) — 7o#(fc — 27r/£) (3) 
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which works only in the scales smaller than the heal- 
ing length £ = 1/f^/g to remove compressible short- 
wavelength excitations. The chemical potential de- 
pends on time to conserve the total number of particles 
N = J dx f(x,t) 2 . The Fourier transformation of the 
moving random potential V{x, t) is V(k, t), which moves 
so that it may satisfy the two-point correlation 



Vq exp 



{x-x'f (t-t'f 



2X1 



2T 2 



{V(x,t)V(x',t')) 

(4) 

where Vq, Xq and To are the characteristic strength and 
scales of space and time, respectively. Thus V{k,t) in- 
jectes the energy to the fluid at the scale Xq to set the 
energy-containing range. 

We numerically solve Eq. J5J in a periodic box with 
spatial resolution containing 256 3 grid points. We con- 
sider the case of g = 1 . As numerical parameters, we use 
a spatial resolution Ax = 0.125 and L = 32, where the 
length scale is normalized by the healing length £. The 
parameters of the moving random potential are Vq = 50, 
X a = 4 and T = 6.4 x 10~ 2 . The numerical integral 
is made by the pseudo-spectral method in space and the 
Runge-Kutta-Verner method in time with a time resolu- 
tion At = 1 x 10 -4 . We start from an initial configura- 
tion in which both the density and the phase of the wave 
function are uniform. The system gets to steady turbu- 
lence after i ~ 25 when both the total energy E and the 
incompressible kinetic energy E^. in , given by 



E =— [ da?$*fV 2 + -/ 2 1$, 
N J 1 2 J 1 ' 
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dx [{/V0}f , 



are almost constant in time. Here {• • • } denotes the in- 
compressible part V •{•••} = of the vector fields. In 
the steady state we can obtain the dissipation rate s from 
the time derivative of the incompressible kinetic energy 
dEL n /dt after switching off the moving random poten- 
tial. 

In this steady turbulence, we can clearly define the 
energy-containing range k < 2t:/Xq in which the sys- 



tem takes the energy from the moving random poten- 
tial, the energy-dissipative range k > 27r/£ in which the 
energy are dissipated by 7(fc), and the inertial range 
2tt/Xq < k < 27r/£ between these two ranges. Figure 
^ (a) shows the energy spectrum of the incompressible 
kinetic energy -E^inCO defined as ^ in = / dfc££ in (fc). In 
the inertial range, the numerically obtained spectrum is 
quantitatively consistent with the Kolmogorov law using 
the above obtained e. Figure Q] (b) shows the flux II(fc) 
of the incompressible kinetic energy which is given by 
the energy budget equation of Eq. @ . The energy flux 
Il(fc) is about independent of k in the inertial range and 
close to the energy dissipation rate e. This result sup- 
ports strongly the scenario of QT played by quantized 
vortices; the energy-containing range, the inertial range 
and the energy-dissipating range are, respectively, made 
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FIG. 1: Dependence of the energy spectrum E^ in (k) (a) and 

(5a) 

the energy flux II(fc, t) (b) of the incompressible kinetic energy 
on the wave number k. Solid lines refer to the Kolmogorov 
(5b) l aw Ce 2 ^ 3 k~ 5 ^ 3 in (a) and the energy dissipation rate e in (b), 
respectively. 



by their nucleation, their Richardson cascade and their 
decay to compressible excitations, as believed in eddies 
of CT. 

The detailed studies of this issue are in progress, being 
reported shortly elsewhere. 
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